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Abstract 

Low-energy particle production perpendicular to the event plane in three-jet 
events produced in Z decays in e + e~ annihilation is measured and compared 
to that perpendicular to the event axis in two-jet events. The topology de- 
pendence of the hadron production ratio is found to agree with a leading-order 
QCD prediction. This agreement and especially the need for the presence of 
a destructive interference term gives evidence for the coherent nature of gluon 
radiation. Hadron production in three-jet events is found to be directly propor- 
tional to a single topological scale function of the inter-jet angles. The slope 
of the dependence of the multiplicity with respect to the topological scale was 
measured to be: 

2.211 ± 0.014(stat.) ± 0.053(syst.) 

in good agreement with the expectation given by the colour-factor ratio 
Ca/Cf = 9/4. This result strongly supports the assumption of local parton- 
hadron duality, LPHD, at low hadron momentum. 
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1 Introduction 

Quantum-mechanical interference effects are basic to the gauge theories of fundamental 
interactions such as Quantum Chromodynamics, QCD. The probabilistic Parton Model 
of strong interactions, however, neglecting interference effects still successfully describes 
many processes as an incoherent sum of the individual sub-processes. In fact it has 
proved difficult to show strong-interaction interference effects in high-energy inelastic 
processes like e + e~ annihilation, deep-inelastic scattering or pp interactions. Evidence 
for coherent soft gluon emission in multiparticle production comes from the necessity 
to include angular ordering in the fragmentation models in order to describe the energy 
dependence of hard interactions, from the so called hump backed plateau in the logarithmic 
scaled-momentum spectrum of hadrons due to the suppression of low-energy particle 
production, and from the string-effect in three-jet events in e + e~ annihilation explained 
by destructive interference. For a comprehensive review see [I]. 

Arguments against the conclusiveness of these verifications of coherence effects have 
been raised, however. Incoherent fragmentation models involving a large number of 
parameters allow for a sufficient description of the data at least at fixed centre-of-mass 
energy [2]. The hadronic momentum spectrum and the energy dependence of the peak 
of this distribution may also be accounted for when assuming a non-minimal phase-space 
structure [3]. The string-effect measurements, except for symmetric topologies [4,1], are 
also influenced by boost effects. 

In this letter evidence for quantum-mechanical interference or colour coherence is pre- 
sented from the comparison of low-energy hadron production perpendicular to the event 
plane of three-jet events and the event axis of two-jet events, respectively. Similarly to 
the string-effect, in this case the presence of destructive interference is expected. Pre- 
suming colour coherence, the measurement is also directly sensitive to the colour of the 
underlying partons, q, q and g. The ratio of the colour charges or the colour-factors for 
gluons and quarks, Ca/Cf, can be verified directly from the dependence of the hadron 
production on the three-jet topology. 

This letter is organised as follows. In Section 2 a basic description of the theoretical 
formulae and the ideas underlying this measurement are given. Section 3 gives an overview 
of the data, selections and corrections applied and describes the measurement. In Section 
4 results on the hadron multiplicity observed in cones oriented perpendicular to the event 
plane of three-jet events and the event axis of two-jet events, respectively, are presented 
and conclusions are given. 

2 Motivation and Theoretical Basis 
of the Measurement 

The hadronisation process is believed to commence with the radiation of soft gluons 
from the initial hard partons. Due to the higher colour charge of gluons compared to 
quarks, soft gluon radiation and consequently hadron production is expected to be in- 
creased in the high energy limit by the gluon to quark colour-factor ratio Ca/Cf= 9/4. 
This simple expectation is perturbed for hadrons sharing a large fraction of the energy 
of the underlying parton. Leading-particle effects become evident and, as such particles 
are created last in time [5], effects due to energy conservation are important. The latter 
influence is stronger for gluon jets due to the more copious particle production. 
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When a radiated hard gluon stays close to the radiating quark, additional soft gluons 
may be unable to resolve the individual colour charges of the close-by partons. In this 
limit only coherent radiation off the parton ensemble, that is from the colour charge of the 
initial quark occurs. This fundamental effect also leads to the angular-ordering property 
of soft gluon emission. 

Coherence effects should be best visible for low-energy hadrons emitted perpendicu- 
larly to the radiating partons as these particles cannot be assigned to a specific jet and 
have a poor resolution power for close-by colour charges. 

The cross-section for soft-gluon radiation in (qqg) three-jet events has been calculated 
in [6]. It varies strongly as a function of the three-jet topology. In contrast to many other 
event characteristics, the soft radiation out of the qqg plane is dominated in the perturba- 
tive calculation by the leading-order contribution which takes a particularly simple form 



with the so called antenna function 

Q, . 

% j = 1 — cos % = 2 sin 2 

Here do~% is the cross-section for soft gluon emission perpendicular to the qqg plane, do2 
the corresponding cross-section for emission perpendicular to the two-jet event (qq) axis, 
and Oij denotes the angle between the partons i and j. The complete topology dependence 
can be factorised in the topological term r t which, multiplied by the energy, obviously 
plays the role of the evolution scale of the low energy hadron multiplicity. 

Equation 1 interpolates between the two cases in which qqg events appear as two-jet 
events. For gluon radiation parallel to the quark (or antiquark) r t = 4/9 and da 3 = da 2 , 
as soft particles cannot resolve the individual colour charges of the quark and the gluon. 
In the case where the gluon recoils with respect to the quark antiquark pair r t — 1. This 
case is similar to an initial two-gluon state forming a colour singlet and dos = CAlCp-da2- 
The term inversely proportional to the square of the number of colours, N c = 3, in r t is due 
to destructive interference and provides an additional suppression of particle production 
on top of the boost effects caused by the string-like topology in general three-jet events. 

The interest of the measurement presented in this paper is: 

• the experimental test of the predicted topology dependence manifesting the coherent 
nature of hadron production, 

• the verification of the 1/iV 2 term as direct evidence for the negative interference 
term beyond the probabilistic parton cascade and 

• the measurement of the slope of the homogenous straight-line given by Equation 1 
and thereby the verification of the colour-factor ratio Ca/Cf- 

In order to verify the above predictions experimentally, charged-particle production 
has been studied in cones with $ con c — 30° half opening angle situated perpendicular 
to both sides of the three-jet plane or the two-jet event axis. For two-jet events the 
azimuthal orientation of the cones is taken randomly. In order to distinguish events with 
two, three and four or more jets the angular ordered Durham (aoD) algorithm is used 
with a jet-resolution parameter y cut = 0.015. The aoD algorithm corresponds to the 
Cambridge algorithm without soft freezing [8]. It has been chosen as it is expected to 
lead to optimal reconstruction of jet directions but to avoid complications due to soft 



[7]: 




F 




(1) 
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freezing. The choice of a fixed jet-resolution is necessary as the leading-order prediction 
of Equation 1, applies to two- and three-jet events. Consequently, events with less than 
two or more than three jets are explicitely excluded from the analysis. 

These experimental choices, made in order to approximate the theoretical prediction, 
imply systematic uncertainties in the data-to-theory comparison. The possible systematic 
variations of the data were assessed by varying the cone opening angle in the range 
20° < $ cone < 40° and the jet-resolution in the range 0.01 < y cut < 0.02. Moreover the 
Cambridge [8] and Durham [9] algorithms were used alternatively. The uncertainty of 
the leading-order prediction, Equation 1, is not specified in [7] and, consequently, is not 
included in the uncertainties of the derived quantities quoted below. 

In order to compare the prediction, Equation 1, to the data directly, the prediction 
is evaluated for each of the possible assignments of the gluon to the three jets. Finally, 
the three results are averaged with each assignment being weighted according to the 
respective three-jet matrix element. The topological term r t in principle varies between 
4/9 ~ 0.44 and 1. In this analysis the lower limit is rj mn ~ 0.5, as the cut in the jet- 
resolution implies a minimum angle between the gluon and a quark. As gluon-jets stay 
unidentified and gluon radiation at low angles dominates, the maximum weighted average 
value of r t , r™ ax ~ 0.71 is obtained for fully-symmetric three-jet events. 



3 Data and Data Analysis 

This analysis is based on hadronic Z decays collected in the years 1992 to 1995 with 
the Delphi detector at the Lep e + e~ collider at CERN. The Delphi detector was a 
hermetic collider detector with a solenoidal magnetic field, extensive tracking facilities 
including a micro-vertex detector, electromagnetic and hadronic calorimetry as well as 
strong particle identification capabilities. The detector and its performance are described 
in detail elsewhere [10,11]. 

In order to select well-measured particles originating from the interaction point, the 
cuts shown in Table 1 were applied to the measured tracks and electromagnetic or 
hadronic calorimeter clusters. Here p and E denote the particle's momentum and en- 
ergy, -^poiar denotes the polar angle with respect to the beam, e k is the distance of closest 
approach to the interaction point in the plane perpendicular to [xy] or along [z] the beam, 
respectively, L t rack is the measured track length. -Ehpc (-^emf) denotes the energy of a 
cluster as measured with the barrel (forward) electromagnetic calorimeter, respectively, 
and -Ehac the cluster energy measured by the hadronic calorimeter. 

The general event cuts shown in Table 2 select hadronic decays of the Z and suppress 
background from leptonic Z decays, 77 interactions or beam-gas interactions to a neg- 
ligible level. Further reduction of background is achieved by the jet-selection cuts given 
also in Table 2. The cut variables are the visible charged energy, ^charged and E^^^ ', 
observed in the event or in one event hemisphere, respectively. Event hemispheres are 
defined by the plane perpendicular to the sphericity axis. The polar angle of this axis 
with respect to the beam is "^sphericity and A^ harged is the observed charged multiplicity. 
Events are discarded, if they contain charged particles with momenta above the kinematic 
limit. 

Events are then clustered into jets using the aoD (or, alternatively, the Durham or 
Cambridge) algorithm. The three-jet event quality requirements also shown in Table 
2 assure well-measured jets. Here 6i denotes the angle between the two jets opposite 
to jet i, i?j et the polar angle of a jet and E^ibie/jet the total visible energy per jet. 
The sequence number % of a jet is given by ordering the jets with decreasing jet energy 
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v^ria hip 


cut 


V 


^ u.4 vjev 


,9 

u polar 


ZX) — 1DU 




^ o cm 


e z 


< 10 cm 


-^track 


> 30 cm 


Ap/p 


< 100% 


-^HPC 


0.5 GeV - 50 GeV 


-Eemf 


0.5 GeV - 50 GeV 


-Ehac 


1 GeV - 50 GeV 



variable 


cut 


general events 


rnhcmisph. 
charged 


> 0.03 • ^/s 


ptotai 
charged 


> 0.12- ^ 


-^charged 


> 5 


^sphericity 


30° - 150° 


Pmax 


45 GeV 


three-jet events 


£i=l °i 


> 355° 


-Evisibic/jet 


> 5 GeV 


Charged /jet 


> 2 




30° - 150° 



Table 1: Selection cuts applied to 
charged-particle tracks and to calorimeter 

clusters Table 2: Selection cuts applied to gen- 

eral events and to three-jet events. 

(Ei > E2 > E3). These exact jet energies were calculated from the inter-jet angles 
assuming massless kinematics (see e.g. [12]). The quality of the selected two-jet events 
is guaranteed already by the event cuts. 

Besides general three-jet topologies, mirror-symmetric three-jet topologies (defined by 
$3 — $2 < 5°) are considered in later cross-checks. Symmetric three-jet topologies form 
a class of especially simple and clean three-jet topologies which we used extensively in 
previous studies [4,12,13]. Therefore, events with symmetric topologies are considered 
separately in order to compare the results obtained with this clean, but statistically 
limited class of events with the results obtained from events with more general topologies. 
Note, that the symmetric-event sample is almost completely included in the general 
sample. However, it allows the angular range to be extended to lower values of Q\. 

Initially 3,150,000 events measured by Delphi enter into the analysis. In total 
1,797,429 Z events fulfill the above event cuts, of which 1,031,080 are accepted as two-jet 
events, 309,227 as general and 53,344 as symmetric three-jet events. 

The charged-particle multiplicity in cones perpendicular to the event axis of a two-jet 
event or the event plane of a three-jet event has then been determined. Besides the cuts 
mentioned above it was required that the axis of the cones has a polar angle of at least 
30° with respect to the beam in order to assure a sufficient acceptance for particle tracks. 

The multiplicity has been corrected for the applied cuts and the limited detector 
acceptance and resolution by a multiplicative correction factor. Using simulated events 
this has been calculated as the ratio of the multiplicities of generated events to accepted 
events at the detector level. The events were generated with the Jetset 7.3 parton- 
shower model [14] as tuned to Delphi data [15]. The influence of the magnetic field, 
detector material, signal generation and digitisation was simulated. The simulated data 
were then treated like the measured events. 

It has already been shown [2,15], that the Monte Carlo model underestimates the 
production of particles at large angles with respect to the primary partons. Therefore, 
the hadron multiplicity in cones perpendicular to the two-jet event axis or the three-jet 
event plane has been reweighted with a small overall multiplicative correction. It was 
deduced from a comparison of data and simulation for two-jet events. After correction 
the model well describes the data and its dependences on the three-jet topology. 

In particular, due to the 400 MeV momentum cut on the charged particles the correc- 
tion factors obtained for the charged multiplicity in the cones are comparably big (~ 2), 
but vary only slightly with event topology. Moreover, in the important ratio of the mul- 
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tiplicity of three- and two-jet events the correction factors cancel to a large extent. The 
correction factor for the ratio N 3 /N 2 is only ~ 0.9. The stability of the acceptance correc- 
tion has been tested by varying the selections imposed on charged-particle tracks shown in 
Table 1. In particular the lower momentum cut has been changed by ±100 MeV. No sig- 
nificant change in the results was observed. Conservatively a 10% relative uncertainty of 
the respective acceptance correction is assumed to propagate to the quoted multiplicities 
or multiplicity ratios. As this uncertainty is expected to mainly change the normalisation 
of the multiplicities it is not included in the overall uncertainties shown in the figures. It 
is, however, considered in the results obtained from global fits of the data. 

4 Results 

Experimental results on the average charged-hadron multiplicities in cones are pre- 
sented for cones perpendicular to the event axis in two-jet events and cones perpendicular 
to the event plane in three-jet events. Moreover the ratio of the cone multiplicities in 
three- and two-jet events is given as well as the momentum spectra in cones perpendicular 
to the plane of three-jet events. 

The charged-hadron multiplicity in cones perpendicular to the event axis in two-jet 
events, N 2 , is given in Table 3 for three different cone opening angles and three different 
values of the jet-resolution parameter y^. The multiplicity shows the expected slight 
increase with y cut and, due to the increase of phase-space, an approximately quadratic 
increase with cone opening angle. 

The corresponding three-jet multiplicity, N 3 , is shown in Figure 1 as a function of 
the inter-jet angles 9 2 and 9 3 for general three-jet topologies and as a function of 9 3 for 
symmetric ones. The cone opening angle is "$ C onc — 30° and = 0.015. For the general 
three-jet events the bin width in 9 2 is: 2° for 98° < 9 2 < 108°, 3° for 108° < 9 2 < 111°, 
4° for 111° < 9 2 < 123° and 5° for 123° < 9 2 < 148°. The corresponding binning in 9 3 
is: 10° for 120° < 9 3 < 140° and 5° for 140° < 9 3 < 165°. For symmetric topologies the 
range in 9\ is fixed by the highest accessible angle of 120° and the bin width of 5° imposed 
by the condition 9 3 — 9 2 < 5°. The lower limit of 9\ ~20° depends on the applied value 
of y cut . Numerical values of the results are given in [16]. The inner error bars represent 
the statistical error, the outer bars also include the systematic uncertainty relevant for 
the comparison to the theoretical expectation of Equation 1, added in quadrature. As 
only the ratio of multiplicities in three- and two-jet events enters in the comparison to 
the theoretical expectation (see Equation 1), the relative systematic uncertainty of the 
multiplicity has been taken to be identical to the relative error of the three- to two-jet 



u cone 


2/cut 


0.01 


0.015 


0.02 


20° 


0.231 ±0.001 


0.245 ±0.001 


0.255 ±0.001 


30° 


0.537 ±0.001 


0.570 ±0.001 


0.593 ±0.001 


40° 


1.007 ±0.002 


1.067 ±0.002 


1.111 ±0.002 



Table 3: Charged-hadron multiplicity in cones perpendicular to the event axis in two-jet 
events, N 2 , as a function of the cone opening angle and the jet-resolution for the aoD 
algorithm. Errors are statistical. 
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multiplicity ratio. It has been determined as the r.m.s. of the results obtained when 
varying the cone size, jet-resolution and jet algorithm as described in Section 2. 

The theoretical expectation, shown as a full line in Figure 1, describes the data well 
for all topologies except for large 6% for symmetric three-jet topologies. This region 
corresponds to small 6\, i.e. close-by low-energy jets. In this case the event plane is ill 
determined and data and theory cannot be compared reliably. 

The prediction increases (by ~ 7%, see dashed line in Figure 1) if the interference term 
of the theoretical prediction oc 1/N% is omitted. In this case the prediction is significantly 
above the data. 

In order to quantify this result the interference term in Equation 1 has been multiplied 
by an amplitude factor k and then fitted to the ratio N 3 /N 2 of the multiplicities as given 
in Figure 1 and Table 3, respectively. The fit yields: 

k = 1.37 ± 0.05 (stat.) ± 0.33 (syst.) X?/Ndf = 1.2 

for general three-jet topologies and 

k = 1.23 ± 0.13 (stat.) ± 0.32 (syst.) X*/Ndf = 2.1 

for symmetric three-jet events with 9\ > 48°. The systematic uncertainty of k was taken 
to be the r.m.s. of the results obtained when varying the cone size, jet-resolution and 
jet algorithm as described in Section 2. Moreover a 10% uncertainty on the acceptance 
correction applied for the multiplicity ratio has been assumed. The x 2 has been calcu- 
lated using statistical errors only. Therefore, a significant destructive interference term 
is observed in the data-to-theory comparison. The size of the term is consistent within 
error with expectation. Moreover it has been verified that the topology dependence of 
the term is consistent with that present in the data. This result represents the 

first direct verification of a colour suppressed negative antenna contribution not included 
in the usual probabilistic parton cascade calculation. Note, that it is based on the di- 
rect comparison of data to an absolute prediction. In particular no construction of an 
unphysical incoherent model is needed for the comparison. 

Figure 2 a) shows the multiplicity ratio N3/N2 as function of r t for different 63. The 
data for different values of #3 agree well as is already expected from the good description 
of the data by the expectation, Equation 1, shown in Figure 1. The exclusive dependence 
on r t as determined from the general-event sample is shown in Figure 2 b) compared 
to the expectation, i.e. to the homogeneous straight-line N 3 /N 2 = Ca/Cf • r t indicated 
dashed. The deviation of the points at low values of r t is again due to events with close-by 
jets. A homogeneous straight line fit to the data for r t > 0.525 yields for the slope: 

2.211 ± 0.014(stat.) ± 0.053(syst.) X?/Ndf = 1.3 

in astonishingly good agreement with the leading-order QCD expectation Ca/Cf= 2.25. 
The systematic uncertainty was again determined as for k in the interference term. The 
fit result is indicated by the full line in Figure 2 b). This result for the first time quanti- 
tatively verifies the colour-factor ratio Ca/Cf from hadronic multiplicities and a leading- 
order prediction. This simple prediction applies, because higher order contributions prove 
to be unimportant for soft particles emitted at large angles [7]. This is different from the 
more global case which depends on all the radiation inside jets. Moreover, these particles 
are least affected by energy- conservation and leading-particle effects. 

For completeness it should be mentioned that the data do not allow the colour-factor 
ratio and the coherence term to be simultaneously verified due to the obvious high cor- 
relation of both terms. If such a fit is attempted the data can as well be described with 
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Figure 1: Multiplicity in cones of 30° opening angle perpendicular to the three-jet event 
plane, N 3 , as a function of the opening angles 2 ,0 3 . The inner error bars are statistical, 
the outer also include systematic uncertainties (see text). The full line is the expectation 
deduced from Equation 1 using N 2 from Table 3. For the dashed line the interference 
term in Equation 1 (oc l/iV c 2 ) has been omitted. 
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r t 

Figure 2: Multiplicity ratio N3/N2 in cones of 30° opening angle as function of r t : a) for 
different angles #3, errors are statistical; b) averaged over # 3 , inner errors are statistical, 
outer errors include systematic uncertainties (see text). The dashed line is the expectation 
Equation 1, the full line is a fit to the data in the indicated r t - interval. 

Ca/Cf^ 2 and N c — > 00. This result is also expected from gauge theories with a large 
number of colours. It is clear, however, that this solution does not apply as Nc = 3 is well 
known experimentally, e.g. from the ratio, R, of the hadronic to the leptonic cross-section 
at the Z pole. 

The dependence of the multiplicity on the scale r t only allows inclusive distributions of 
the produced hadrons to be studied. In Figure 3 a) the differential momentum distribution 
is shown for different values of r t . The momentum distributions are observed to scale 
approximately with r t . This is most clearly seen in Figure 3 b) where the ratio of the 
momentum spectra to the average spectrum is shown. Moreover, these ratios have been 
scaled by (r t ) / r t such that a unit value for the overall multiplicity ratio is expected. For 
small momenta (p < 1 GeV) the data are consistent with unity, for higher momenta less 
(more) particles are produced for small (high) scales r t , respectively. The constancy of the 
ratio is expected if the production of hadrons is directly proportional to gluon radiation, 
thus the observed result strongly supports the conjecture of local parton-hadron duality, 
LPHD [17], for low-energy hadrons. The deviation from unity for higher momenta may 
be understood as being due to the increase of the hadronic phase-space for high r t . 
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Figure 3: a) Differential momentum distribution of hadrons in cones of 30° opening 
angle perpendicular to the three-jet event plane for different values of r t (aoD algorithm, 
Ucut = 0.015). b) Ratio of the momentum distributions to their average reweighted using 
Equation 1 and the individual r t values. 



In summary the multiplicity in cones perpendicular to the event plane of three-jet 
events produced in e + e~ annihilation has been measured for the first time and com- 
pared to the corresponding multiplicity in two-jet events. The multiplicity ratio is well 
described by a leading-order QCD prediction also including a destructive interference 
term. The presence of destructive interference and thereby of the coherent nature of soft 
hadron production was verified from the comparison of the data to the absolute QCD 
prediction. The slope of the hadron multiplicity with the topological variable r t directly 
represents the colour-factor ratio Ca/Cf- This measurement for the first time is able to 
verify the colour-factor ratio using a hadron multiplicity measurement and a leading-order 
QCD prediction. This possibility is a consequence of the favorable perturbative situation 
studied and the focus on low-energy large-angle particles. A simultaneous study of the 
momentum and topology dependence of the hadron multiplicity provides evidence for the 
LPHD conjecture. 
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